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Using general arguments of an optimization taking place between the pair wave function and the 
repulsive part of the electron-electron interaction, we analyze the superconducting gap in materials 
with multiple Fermi-surface (FS) pockets, with exemplary application to two proto-type ferrop- 
nictide setups. On the basis of functional renormalization group (FRG) calculations for a wide 
parameter span of the bare interactions, we show that the symmetry of the gap and the nodal 
versus nodeless behavior is driven by this optimization requirement. 
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The discovery of superconductivity (SC) in new com- 
pounds such as the iron pnictides lias opened up a new 
avenue for studying the fundamental question "what is 
universal and what is material-dependent" concerning 
the mechanism of high-Tc SC in a wider class of materi- 
als other than, but also including, the cuprates. There, 
after more than two decades of intense research, more 
and more theoretical as well as experimental studies sup- 
port a scenario where the general nature of the rf-wave 
SC as well as other salient features are accounted for 
by an electronic pairing mechanism extracted from a 
one-orbital Hubbard model [TJ [2] with the addition that 
the material-dependence is embedded in the multi-orbital 
(e.g. 3-band) extensions [5]. 

In other SC compounds such as the pnictides, however, 
the picture seems more complicated: Here, at the outset, 
multi-band SC appears with gaps possibly displaying dif- 
ferent symmetries such as extended (sign-reversing) s- 
competing with d-wave and with nodal or also nodeless 
behavior on the disconnected Fermi surface (FS) sheets 
[iHTSl \TT[ [55]. Accordingly, even the simplest multi- 
band Hamiltonian with only on-site interactions contains 
four possibly relevant terms, the intra-orbital and inter- 
orbital repulsion as well as the Hund's-rule coupling and 
pair hopping. Searching for SC pairing, these interac- 
tions have to be augmented with the orbital dependence 
of the FS pockets, since the interactions become ma- 
trices formed by local orbitals which have a dominant 
orbital "weight" at the FS pockets. We hence investi- 
gate whether this intricate interplay of multi-orbital band 
structure, FS topology and interactions still allows for 
insights into a more universal than material-dependent 
understanding of SC in these systems. 

In this Letter, we describe an "optimization princi- 
ple" which can help providing such a more universal pic- 
ture. The main point of our work is to show that the SC 
state, its gap, and, in particular, its anisotropy in mo- 
mentum space is determined by an optimization which 
determines and optimizes the interplay between the at- 



tractive interaction in the SC-channel and the Coulomb 
repulsion. This optimization problem, as discussed be- 
low, is unavoidable in a multi-band SC situation: for the 
pnictides, it appears because of a frustration in the s±- 
channel, when more than two FS-pockets are involved in 
setting up the pairing interaction. 

Already from the BCS gap equation, one can see that 
Coulomb repulsion at a finite momentum transfer can in- 
duce pairing only when the wave vector of such an inter- 
action connects regions on one FS (in the cuprate case), 
or regions on different FSs (in the pnictide case), which 
have opposite signs of the SC order parameter. This cor- 
responds to putting the electron pairs in an anisotropic 
wave function such as d-wave in the high-Tc cuprates, 
or the sign- reversing s-wave {s± ) in the pnictides, where 
in the latter case the wave vector (7r,0) in the unfolded 
Brillouin Zone connects hole (h) and electron (e) FS- 
pockets with a sign-changing s± gap [U [8] . Early studies 
based on either RPA spin-fluctuations (SF) scenarios [4] 
or on Renormalization-Group [RG] studies of just one- 
hole and one-electron FS have reported a momentum- 
independent s± gap. At first glance, this similarity of 
the gap function obtained by so dissimilar approaches 
as RPA and FRG may appear surprising. Indeed, the 
repulsive part of the Coulomb interaction is treated dif- 
ferently which leads to differing results for the general 
multi-pocket case |17) . 

The interesting setup for the optimization principle 
concerns a multi-pocket situation - as generally appear- 
ing in the ferro pnictides - where more than two pockets 
create crucial pairing interactions. In order to illustrate 
the principle at work for such a scenario, we investigate 
a 4-pocket and a 5-pocket Fermi surface (pFS) topology 
originating from a 5-band model (Fig. [T]), and discuss the 
superconducting order parameter from that perspective. 

The principal physical content of this optimization sce- 
nario can already be observed for the 4pFS and 5pFS 
scenarios in Figs. [T^ and b. Let us try to understand 
the systems from the unfrustrated s± limit where the 
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FIG. 1. (Color online). Scan of the intra-orbital interaction Ui for both the 4pFS (Fig. [T^) and 5pFS (Fig. [T]d) situations. 
Colors on the FS-pockets give dominant orbital weights, with dashed segments signalling several leading contributions. The 
increase of Ui leads to an increase of the e-pocket gap for the 5pFS scenario, while it favors the nodal scenario for 4pFS. 



r o X pair scattering between h-pockets at r(0,0) and 
3-pockets at X{Tr, 0) is minimized. Here, a green arrow 
foT X X scattering (4pFS, lower FS display in Fig.[I^) 
and a magenta arrow (5pFS) for X ^ M indicate addi- 
tional interactions (dependent on the dominant orbital 
weights on the FSs), The green interaction in Fig. [l^ 
frustrates the previous pure s± limit. The system then 
strikes a compromise - and this is the essence of the op- 
timization principle - by enhancing the anisotropy of the 
gap function (denoted by /^^(k) in Fig. [ij on the e- 
pockets at X (FS positions 32 to 64 in Fig. [T]), even- 
tually reaching even a nodal situation for larger inter- 
actions. The dominant part of the X -^r^ X interaction 
acts so as to push the peaks of the e-gap function fur- 
ther up, while the dominant part of the (F o X) in- 
teraction (magenta arrow in Fig. [l^) tries to push the 
e-gap valleys down. Thus, a transparent understanding 
of the anisotropics and the nodeless versus nodal behav- 
ior emerges: the multi-band SC adjusts the momentum 
dependence of the gap, i.e. its anisotropy, so as to mini- 
mize the effect of the Coulomb repulsion 'TS]. 

In more mathematical terms, this optimization is re- 
flected in Eq. [s] below for the dominant Cooper-channel 
eigenvalue cf^%A.) taking the largest negative value: 

(A) = (/SC(k)y|C(k,_k,p)/SC(p)*) (1) 

Here, as detailed in Eq. [s] V^^ denotes the pairing func- 
tion, where A is the RG-flow parameter and /^'-'(k) the 
SC (gap) form factor associated with it. (...) denotes the 
inner product and involves the k- and p-points on all 4 
(or 5) FS-pockets (Fig.[l]). We have 
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and its largest negative value is determined via an opti- 
mization taking place between all pockets I and to. This 



is a frustration problem as not all minimization condi- 
tions can be fulfilled at the same time. 

In conjunction with the underlying FS topology ob- 
tained from LDA-type calculations fS], our FRG-studies 
also serve to answer which of the four intrasite interac- 
tions in the starting Hamiltonian are playing a leading 
role. To this extent, we report on an extensive parameter- 
sweep study the results of which are fully in line with the 
optimization argument. Placing the value of these inter- 
actions finally around the spread of values obtained in 
recent ab-initio DFT work [5D], we find the intra-orbital 
interaction Ui to take on the pivotal role. The 4pFS and 
5pFS scenarios can be cast into a 5-band model, with 
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Here c's stand for electron annihilation operators, a, b for 
the d-orbitals, and s denote the spin indices. 

In the many-body part the intra- and inter-orbital in- 
teractions Ui and U2, as well as the Hund's coupling Jh 
and the pair hopping Jpair enter, i.e. 
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where 



Hi 



denote density operators at site i of spin s in 



orbital a. Typical interaction settings are dominated by 
intra-orbital coupling, Ui > U2 > Jh Jpair, and can be 
obtained from constrained RPA calculations [20j . 

Details of the electronic structure such as the FS topol- 
ogy and, more specifically, the presence of the (7r,7r)- 
pocket are crucial for the SC state [iHT5]. Using FRG on 
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the above 5-band (Fe d-orbital) model of the Fe-based 
SC with orbital interactions as in Eq. |3] and |4j we have 
recently found that the gap on the e-pockets can un- 
dergo a nodal transition if the h-pocket at (tt, tt) is ab- 
sent [T5] . Similar conclusions have been reached by Maiti 
and Chubukov using a parquet RG analysis [21] ■ On the 
basis of RPA calculations, Kuroki et al. [S] using a sim- 
ilar Hamiltonian as in Eqs. [3] and |4] have already argued 
that the (tt, 7r)-pocket is sensitive to the lattice structure 
("pnictogen height") and crucial for the gap structure. 
Kemper et al. [TB], within again a 5-orbital RPA theory, 
further substantiated these conclusions. In particular, 
they demonstrated the sensitive dependence of the SC 
state to aspects of the electronic structure such as the 
FS topology and the FS orbital weights. 

In the FRG pflU [HI [23], one starts from the bare 
many-body interaction (Eq. (2)) in the Hamiltonian and 
the pairing is dynamically generated by systematically 
integrating out the high-energy degrees of freedom in- 
cluding important fluctuations (magnetic, SC, screening, 
vertex corrections) on equal footing. This differs from 
the RPA which takes right from the outset a magneti- 
cally driven SF-type of pairing interaction. 

For a given instability characterized by some order pa- 
rameter Ok, the 4-point function (4PF) VA(ki, k2, ks, k4) 
in the particular ordering channel can be written in short- 
hand notation as p VA(k, p)[Ojj. Op] [H]. Accordingly, 
the 4PF VA(k, — k, p, — p)) in the Cooper channel can be 
decomposed into different eigenmode contributions [9lll5j 
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where i is a symmetry decomposition index, and the lead- 
ing instability of that channel corresponds to an eigen- 
value c^'-^(A) first diverging under the flow of A. /^*-''*(k) 
is the SC form factor of pairing mode i which tells us 
about the SC pairing symmetry and hence gap struc- 
ture associated with it. In FRG, from the final Cooper 
channel 4PFs, this quantity is computed along the dis- 
cretized Fermi surfaces (as shown in Fig. [T]) , and the 
leading SC instabilities are plotted in Figs, [^o [3] If not 
stated differently, the interaction parameters not speci- 
fied in the plots are kept fixed at the representative setup 
Ui = 3.5ey, U2 = 2.0eV, Jh = Jpair = 0.7eV. 

We first investigate the behavior upon the variation of 
the intraorbital interaction scale Ui (Fig. [T]). For small 
values of Ui{Ui = U2 = 2eV) in Fig. the s±-sign 
change is induced by the "optimization principle" be- 
tween the h-pockets around the F-point and the e-pockets 
around the X points. By increasing Ui, the 4pFS sys- 
tem develops a pronounced gap anisotropy at the electron 
pockets, which eventually leads to gap nodes at C/i ~ 3eV 
(Fig. [1^). This behavior is due to an enhanced Ui re- 
pulsion within the d^^-Y^ orbitals, which amplifies the 
pair scattering within the (ix2-Y2-dominated parts of one 
electron pocket to the other (see green arrow in Fig. [l^ 
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FIG. 2. (Color online). Scan of U2- Upon increase, more 
and more relevant scattering contributions from the 3rd hole 
pocket (M) in the 5pFS scenario (row (b)), to parts of 
the electron pockets of different orbital content increase the 
anisotropy. In the 4pFS scenario (row (a)), the anisotropy 
decreases due to a weakening of the previously dominant scat- 
tering between the peaks and valleys of the e-pockets. 



for Ui = 5eV). The above mentioned optimization re- 
quirement between a repulsive interaction and the pair 
wave function then favors an increased gap anisotropy 
between the peak and the valley of the SC gap on the 
e-pockets, which eventually yields a sign change. 

In the 5pFS scenario (Fig. [T]d), the additional AI- 
pocket, which exclusively carries dx'^^y^ orbital weights, 
also generates pair scatterings due to Ui to the dx^-y^- 
dominated parts of the electron pockets (see magenta ar- 
rows in Fig.[T|D). This pushes down the peaks (tips) of the 
SC-gap on the electron pockets, on the basis of the same 
orthogonality argument (magenta arrow in Fig. [l]3). 

As the inter-orbital interaction U2 is increased (Fig. [2| , 
the main change is that the significance of the dx^-Y^ 
scatterings driven by Ui is slightly lowered. As the in- 
terorbital scattering phase space becomes important, the 
orbital distribution along the pockets determines more 
and more the behavior. For the 5pFS scenario (row (b) 
in Fig. [2| , this gives an increased e-pocket anisotropy and 
a smaller e-pocket gap. For the 4pFS scenario, the nodal 
propensity is significantly reduced, as the previously de- 
cisive scattering between the peaks and the valleys of the 
electron pockets becomes less relevant. 

As Jh is increased, similar to Wang et al. [13], we ob- 
serve that the anisotropy on the e-pockets is enhanced, 
again because of the general orthogonality requirement 
between the repulsive interaction and the SC pair state, 
which applies both to 5pFS and 4pFS scenarios. Within 
a reasonable parameter range up to ^ leV, the modifi- 
cation of the SC form factor is comparably small. 

On the bare level, Jpair gives a positive semidefinite 
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FIG. 3. (Color online). Scan of Jpair. The increase of pair- 
hopping helps to drive a gap for the 4pFS. 



contribution to the Cooper channel. As a consequence, it 
increases the absolute Cooper channel scale and leads to 
an increase of the A-divergence scale of the SC instability 
(not shown here). While the SC form factor for the 5pFS 
scenario remains essentially unchanged, we observe some 
decrease of nodal propensity in the 4pFS scenario (see 
Fig. [3]). There, the Jpair scattering contribution between 
the r hole pockets and the tips of the electron pockets 
play a similar role as the scattering contribution of the 
5th hole pockets in the 5pFS scenario. However, this 
scattering takes place now between dxz — > dx^-y^ and 
dyz ^ dx^-Y^: orbitals (red arrows in Fig. [s]). 

In summary, we have demonstrated the usefulness of 
the optimization principle to provide a more univer- 
sal characterization of gap anisotropies in multi-band 
SC. The optimization minimizes the Coulomb repulsion, 
which is needed because of frustration (such as occurring 
in the s± -channel in our example), when more than two 
FS-sheets are contributing to the pairing interaction. Us- 
ing FRG for two generic FS topology setups of the pnic- 
tides, we show that the optimization principle relates the 
gap anisotropies and their nodeless versus nodal behav- 
ior in a rather transparent way to the multi-orbital band 
structure, FS topology, and interactions. 



ACKNOWLEDGEMENTS 

Useful discussions are gratefully acknowledged with D. 
J. Scalapino, C. Honerkamp and B. A. Bernevig. The 
work was supported by DFG-SPP 1458/1 (CP, RT and 
WH) and by the Bavarian KONWIHR Program. RT is 
supported by the Humboldt Foundation. 



[1] P. W. Anderson, Science 235, 1196 (1987) 
[2] See, for example, D. J. Scalapino, Chapt. 13 in "Hand- 
book of High Temperature Superconductors" editors J. 

R. Schrieffer and J. S. Brooks, Springer (2007), [cond- 

Q mat/0610710 and refs. therein. 
[3] For a recent review, see also "Towards a Better Under- 
standing of Superconductivity at High Transition Tem- 



[5 

[6: 

[7 
[8 
[9 
[10 

[11 
[12 
[13 
[14 
[15 
[16 

[17 



[18] 



[19] 

[20] 

[21] 
[22] 

[23] 

[24] 



pereatures" editors R. Hackl and W. Hanke, Eur. Phys. 
J. Special Topics 188, 3 (2010) 

I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008) 

K. Kuroki, H. Usui,S. Onari, R. Arita, and H. Aoki, Phys. 
Rev. B 79, 224511 (2009) 

S. Graser, T. A. Maier, P. J. Hirschfeld, and D. J. 
Scalapino, New Journal of Physics 11, 025016 (2009) and 
Phys. Rev. B 79, 224511 (2009) 

V. Stanev, J. Kang, and Z. Tesanovic, Phys. Rev. B 78, 
184509 (2008) 

A. V. Chubukov, D. V. Efremov, and I. Eremin, Phys. 
Rev. B 78, 134512 (2008) 

F. Wang, H. Zhai, Y. Ran, A. Vishwanath, and D.-H. 

Lee, Phys. Rev. Lett 102, 1047005 (2009) 

C. Piatt, C. Honerkamp, and W. Hanke, New J. Phys. 

II, 055058 (2009) 

R. Thomale, C. Piatt, J. Hu, C. Honerkamp, and B. A. 
Bernevig, Phys. Rev. B 80, 180505 (2009) 
H. Zhai, F. Wang, and D.-H. Lee, Phys. Rev. B 80, 
064517 (2009) 

F. Wang, H. Zhai, and D.-H. Lee, Phys. Rev. B 81, 
184512 (2010) 

A. V. Chubukov, M. G. Vavilov, and A. B. Vorontsov, 
Phys. Rev. B 80, 140515 (2009) 

R. Thomale, C. Piatt, W. Hanke, and B. A. Bernevig, 
arXiv:1001.3599 (2010) 

A. F. Kemper, T. A. Maier, S. Graser, H.-P. Cheng, P. J. 
Hirschfeld, and D. J. Scalapino, New Journal of Physics 
12, 073030 (2010) 

This similarity is partly based on the fact that the 
anisotropic wave function completely avoids the strongly 
repulsive (short-range) part of the Coulomb interaction 
(i.e. C/ X]fc(cfctC-fct) ~ 0)1 ^ single-band Hubbard-U 
model, or, to a large extent for a 2-pocket toy model |24) . 
This is crucial for RPA-treatments, where the bare re- 
pulsion U appears also in the pairing channel. In the 
RG treatments for a 2-pocket situation the correspond- 
ing intra-pocket repulsive part of the Coulomb interac- 
tion is renormalized in the flow eventually even changing 
its sign to become attractive. Yet, the pairing, i. e. the 
corresponding flow, is dominated by the repulsive inter- 
pocket pair scattering with wave vector (tt, 0) connecting 
the e- and h-pocket 8, 10 . Thus, in the RG, the pair- 
wave function is "optimized" , i.e. gives the largest pairing 
eigenvalue (see below), again for an s±-form. 
Some of the arguments in our Letter for the multi-pocket 
scenarios are independent of the method chosen for ex- 
tracting the pairing and similar to arguments found e.g. 
in Ref. |211I22| . However, the "optimization", i.e. balanc- 
ing the attractive s± interaction and the Coulomb repul- 
sion (frustration) is different in RG versus RPA treat- 
ments (see also Ref. [T7| and |21p. 

For a recent experimental work and refs. therein see e.g. 

H. Kim et al., arXiv:1008. 3251^^3 (2010) 

For a recent review see M. Imada and T. Miyake, ch.5.2.6, 
J. Phys. Soc. Jpn. 79, 112001 (2010) 
S. Maiti and A. B. Chubukov, [arXiv:101 0.0984 v2 (2010) 
T. A. Maier, S. Graser, D. J. Scalapino, and P. Hirschfeld, 
Phys Rev. B 79, 224510 (2009) 

C. Honerkamp, M. Salmhofer, N. Furukawa, and T. M. 
Rice, Phys. Rev. B 63, 035109 (2001) 

I. Mazin and J. Schmahan Physica C, 469, 614 (2009) 



